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ABSTRACT 

Chandra observations of large samples of galaxy clusters detected in X-rays by ROSAT provide a new, robust 
determination of the cluster mass functions at low and high redshifts. Statistical and systematic errors are now 
sufficiently small, and the redshift leverage sufficiently large for the mass function evolution to be used as a use- 
ful growth of structure based dark energy probe. In this paper, we present cosmological parameter constraints 
obtained from Chandra observations of 37 clusters with (z) = 0.55 derived from 400 deg^ ROSAT serendipitous 
survey and 49 brightest z « 0.05 clusters detected in the All-Sky Survey. Evolution of the mass function be- 
tween these redshifts requires Ha > with a ~ 5ff significance, and constrains the dark energy equation of state 
parameter to Wq = -1.14 ± 0.21, assuming constant w and flat universe. Cluster information also significantly 
improves constraints when combined with other methods. Fitting our cluster data jointly with the latest super- 
novae, WMAP, and baryonic acoustic oscillations measurements, we obtain Wg = -0.991 ± 0.045 (stat) ±0.039 
(sys), a factor of 1.5 reduction in statistical uncertainties, and nearly a factor of 2 improvement in systematics 
compared to constraints that can be obtained without clusters. The joint analysis of these four datasets puts a 
conservative upper Hmit on the masses of light neutrinos, £ < 0.33 eV at 95% CL. We also present updated 
measurements of Ciuh and as from the low-redshift cluster mass function. 

Subject headings: cosmology: observations, cosmological parameters, dark matter — clusters: general — surveys 



1. DARK ENERGY AND CLUSTER MASS FUNCTION 

Recent accelerated expansion of the Universe detected in the 
Hubble diagram for distant type la supernovae is one of the 
most signi ficant discoveries of the past 10 years ( Perlmutter] 
|et aL|1999||Riess et al. i 1998 1 . The acceleration can be attributed 
to the presence of a significant energy density component with 
negative pressure, hence the phenomenon is commonly re- 
ferred to as Dark Energy. For a recent review of the dark energy 
discovery and related the oretical and observational issues, see 
[Frieman, Turner & Hutere r ( 2008 ) and references therein. Per- 
haps the simplest phenomenological model for dark energy 
is non-zero Einstein's cosmological constant. The supernovae 
data indicated (and other cosmological datasets now generally 
agree) that a cosmological constant term currently dominates 
energy density in the Universe. 

The next big question is whether Dark Energy really is the 
cosmological constant. The properties of dark energy are com- 
monly characterized by its equation of state parameter, w, de- 
fined as p = wp, where p is the dark energy density and 
p is its pressure. A cosmological constant in the context of 
General Relativity corresponds to a non-evolving w = -1. It 
is proposed that departures from the cosmological constant 
model should be sought in the form of observed w being either 
+ -1, or evolving with redshift. Combination of supernovae, 
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cosmic microwave background, and baryonic acoustic os cilla 



tions data curren tly constrain |1 -i- w| < 0.15 at 95% CL iKo 



matsu et al. 2008} . Observational signatures of such deviations 



ot w from - 1 are very small, and hence the measurements are 
prone to systematic errors. For example, variations of w be- 
tween -1 and -0.9 change fluxes of z = 0.75 supernovae in a 
flat universe with Om = 0.25 by only 0.03 magnitudes. There- 
fore, it is crucially important that the dark energy constraints at 
this level of accuracy are obtained from combination of several 
independent techniques. This not only reduces systematics but 
also improves statistical accuracy by Ijreaking degeneracies in 
the cosmological parameter constraints. 

One of the methods that has been Httle used so far is evolu- 
tion in the number density of massive galaxy clusters. Evolu- 
tion of the cluster mass function traces (with exponential mag- 
nification) growth of linear density perturbations. Growth of 
structure and distance-redshift relation are similarly sensitive 
to properties of dark energy, and also are mutually highly com- 
plementary methods (e.g., Linder & Jenkins 20031. Mapping 
between the linear power spectrum and cluster mass func- 
tion relies on the model for nonlinear gravitational collapse. 
This model is now calibrated extensively by N-body simula- 
tions (see §[3]l. The cluster mass function models also use ad- 
ditional assumptions (e.g., that the mass density is dominated 
by cold dark matter in the recent past, and that the fluctua- 
tions have Gaussian distribution). However, corrections due 
to reasonable departures from these assumptions are neghgi- 
ble compared to statistical uncertainties in the current samples 
(we discuss these issues further in §[3|. It is important also 
that the theory of nonhnear collapse is insensitive to the back- 
ground cosmology. For example, the same model accurately 
describes the relation between the linear power spectrum and 
cluster mass function in the Om = 1> ^^a = 0, low-density 
D.M = 0.3, Oa = 0, and "concordant" Om = 0.3, Ha = 0.7 



cosmologies (Jenkins et al. 2001 1. 
Fitting cosmological mo3els to the real cluster mass func- 



2 



VIKHLININ ET AL. 



tion measurements uses not only growth of structure but also 
the distance-redshift information because observed proper- 
ties for objects of the same mass generally depend on the dis- 
tance. Therefore, constraints on w derived from the cluster 
mass function internally make a combination of growth of 
structure and distance based cosmological tests, and thus po- 
tentially can be very accurate a nd competitive with any other 
technique (e.g.,|Albrecht et al. 2006| ). 

Previous attempts to use evolution of the cluster mass func- 
tion as a cosmological probe were limited by small sample sizes 
and either poor proxies for the cluster mass (e.g., the total X- 
ray flux) or inaccurate measurements (e.g., temperatures with 
large uncertainties). Despite these hmitati ons, reasonable con- 
straints coul d still be derived on Om (e g-> Borgani et al. 2001 
Henry|20"04 l. However, constraints on the dark energ y equa 



tion of state from such studies are weak. For example, Henry] 
( 2004| l derived the best-fit w = -0.42, only marginally incon- 
sistent with w = - I, using the temperature function of the Ein- 



stein Medium Sensitivity Survey clusters; Mantz et al. 1 2007 1 
determine w = - 1 .4 ± .55 with a lar ger sample of distant cius- 
ters (MACS survey, see |Ebeling et al . 200 1 ) but using the X-ray 
luminosity as a mass proxy. 

The situation with the cluster mass function data has been 
dramatically improved in the past two years. A large sample 
of sufficiently massive clusters extending to z ~ 0.9 has been 
derived from ROSA T PSPC pointed data covering 400 deg^ 



([Burenin e t al.|2007 Paper 1 hereafter). Distant clusters from 
the 400d sample were then observed with Chandra, providing 
high-quality X-ray data and much more accurate total mass 
indicators. Chandra coverage has also become available for 
a complete sample of low-z clusters originally derived from 
the ROSAT All-Sky Survey. Results from deep Chandra point- 
ings to a number of low-z clusters have significantly improved 
our knowledge of the outer cluster regions and provided a 
much more reUable calibration of the Mtot vs. proxy relations 
than what was possible before. On the theoretical side, im- 
proved numerical simulations resulted in better under stand- 



ing of measurement biases in the X-ray data analysis INagai 
let al |2007| |Rasia et al . '2006' 'Jelte ma et al. |20 07 ) . Even more 
importantly, results from these simulations have been used to 
suggest new, more reliable X-ray proxies for the total mass 
pCravt sov et al. 2 006 ^ 1. We dis cuss all this issues in the previous 



paper ( Vikhlini n et ai.|2008 Paper II hereafter). The cluster 
veoi 



mass functions derived in this paper are reproduced in Fig.fT] 
Overall, these results are an important step forward in proviS 
ing observational foundation for cosmological work with the 
cluster mass functions. 

In this work, we present c osmological constraints from 
the data discussed in [Paper iT The cosmological information 
contained in the cluster mass function data and relevant to 
dark energy constraints can be approximately separated into 
3 quasi-independent components: 

(1) Changes in the comoving number density at a fixed 
mass threshold constrain a combination of the perturbations 
growth factor and relative distances between low and high-z 
samples; this by itself is a dark energy constraint (^. 

(2) The overall normaUzation of the observed mass func- 
tion constrains the ampUtude of linear density perturbations 
at z « 0, usually expressed in terms of the ffs parameter Sta- 
tistical and systematic errors in the Cg measurement are now 
sufficiently small, and the ratio of ffg and the ampUtude of the 
CMB fluctuations power spectrum gives the total growth of 
perturbations between z « 1000 and z = — a second power- 
ful dark energy constraint (§|8.1|l. 



(3) The slope of the mass function measures Om ^ h; this 
by itself is not a dark energy probe but can be used to break 
degeneracies present in other methods. 

Our dark energy constraints were derived for the follow- 
ing cases. Assuming constant w and flat universe, we mea- 
sure Wo = -1.14 ± 0.21 using only cluster data (i.e., evolution 
of the mass function between our two redshift samples) and 
the HST prior on h (§8.21. Combining cluster and WMAP 
data, we obtain Wq = - 1 .08 ± 0.15 but (wq is constrained much 
more tightly for a fixed CIm (§ |8.3| . Finally, adding cluster data 
to the joint supernovae + WMA P -i- BAO constraint, we ob- 
tain Wq = -0.991 ± 0.045 (§I83|, sign ificantly reducing statis- 
tical and especially systematic (§ |8.4| uncertainties compared 
to the case without clusters. A large fraction of the extra con- 
straining power comes from contrasting as, with normaliza- 
tion of the CMB power spectrum; this procedure is sensitive 
to non-zero mass of Ught neutrinos. Allowing for > 0, 
we obtain a new conservative upper limit Y,'^v < 0.33 eV 
(95% CL) while still improving the Wq measurement relat ive to 
the SN-hWMAP-hBAO-only case (wo = -1.02 ± 0.055, §[85). 
Adding clusters also improve s equ ation of state constrains for 
evolving w in fl at universe (§ |9.1| and constant w in non-flat 
universe (§ |9.2| l 

The paper is organized as follows. We start with a short sum- 
mary of cluster data and systematic uncertainties (§[2j, discuss 
issues relevant for computing theoretical mass function mod- 
els (§[3} and describe our fitting procedure (§[5}. We then dis- 
cuss constraints that can be obtained from low-redshift mass 
function only {Cluh in §|5]and Cg in §[6|. We then consider as 
an example constraints from the cluster evolution in non-flat 
ACDM model (i.e., w fixed at -1); CIa > is required with 
~ 5a confidence (§[7|. Constraints on the dark energy equa- 
tion of state are c onsidered in §§[8]-|9] Systematic errors are 
discussed in § |8.4| 

2. SUMMARY OF THE CLUSTER DATA AND 
SYSTEMATIC UNCERTAINTIES 

This work is based on two cluster s amples, o riginally com- 
piled from ROSAT X-ray surveys (see Paper II| for a complete 
description of the sample selection and data analysis). The 
low-redshift sample includes the 49 highest-flux clusters de- 
tected in the All-Sky Survey at Galactic latitudes \b\ > 20° 
and z > 0.025. The effective redshift depth of this sample is 
z < 0.15. The high-redshift sample includes 37 z > 0.35 ob- 
jects detected in the 400d survey, with an additional flux cut 
applied; the redshift depth of this sample is z r* 0.9. All the low 
and high-z clusters were later observed with Chandra, provid- 
ing good statistical precision spatially resolved spectral data 
thus yielding several high-quality Mtot estimators for each ob- 
ject. The combined cluster sample is a unique, uniformly ob- 
served dataset. The volume coverage and effective mass Umits 
at low and high redshifts are similar (see the estimated mass 
functions in Fig.[T]l. 

Because of the sufficiently high quaUty of the Chandra data, 
we employ advanced data analysis techniques going well be- 
yond simple flux estimates and /? -model fits commonly used 
in earlier studies. Cosmological cluster simulations has been 
used to test for the absence of significant obser vational bi- 
ases i n reconstructing the basic cluster parameters ( Nagai et al. 
2007| . Using these simulations, we also tested which of the 
X-ray observables are best proxies fo r the total cluster mass 
I Kravtsov et al.|2006|[Nagai et al.|2007) l and concluded that the 
best three are the average temperature, Tx, measured in the an- 
nulus [0.15-1] rsoo (thus excluding the central region often af- 
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Fig. 1. — Estimated mass functions for our cluster samples computed for the 
Om = 0.25, Cl\ = 0.75, h = 0.72 cosmology. Solid lines show the mass 
function models (weighted with the survey volume as a function of M and z), 
computed for the same cosmology with only the overall normalization, (is , fit- 
ted. The deficit of clusters in the distant sample near Msoo = 3 X 10''* ft"' Mq 
is a marginal statistical fluctuation — we obs erve 4 clusters where 9.5 are ex- 
pected, a 2a deviation (c£ Fig. 17 in jPaper fl ;). 



fected by radiative cooling and sometimes, by AGN activity in 
the central galaxy); the intracluster gas mass integrated within 
rsoo; and the combination of the two, Yx = Tx x Mgas. These 
parameters are low-scatter proxies of the total mass (in par- 
ticular, Yx, and Mgas is only slightly worse). Simulations and 
available data show that the scaUng of these proxies with Mtot, 
including the redshift dependence, is very close to predictions 
of the simple self-similar model. In a sense, even though we 
use advanced numerical simulations which include multiple 
aspects of the cluster physics to test Mtot vs. proxy relations, 
the role of simulations is Umited to providing small corrections 
to predictions of very basic and hence reUable theory. AppU- 
cation of these corrections as well as practical considerations 
fo r deriving Tx, Mgas, and Yx from the real data are discussed 
in|Paper III 

Paper II also presents an observational calibration of the 
Mtot vs. proxy relations using an extremely well-observed sam- 
ple of low-z clusters. This discussion is crucial for understand- 
ing the systematic uncertainties in our cluster mass func tion 
measu rements, and we urge interested readers to consult 



Pa- 



per II Table 4 there gives a summary of the main sources ot 
systematic uncertainties in the derived cluster mass functions. 
They can be separated into three quasi-independent compo- 
nents. First is the uncertainty in calibration of the absolute 
cluster mass scale by Chandra hydrostatic mass estimates in 
a sample of dynamica lly relaxed, well-observed low-z clusters 
( Vikhlinin et al.|2006f ; the level of this uncertainty (9%) is es- 
timated from compar ison of Chandra masses with two recent 
weak lensing studies (Hoekstra 2007 Zhang et al. 2008). Sec- 
ond is uncertainties related to possible departures from stan- 



dard evolution in M,, 



Tx, Mu 



Mgas, and Mtot - ix 



Yx rela- 



tions. This uncertainty (~ 5 - 6% between z = and z = 0.5) 
was estimated from general reliability of numerical models of 
the cluster formation and from the magnitu de of corr ections 
that had to be applied to the data (see §4 in [Paper iTlfor de- 



tails). The last major source of uncertainty is evolution in the 
Lx - Mtot relation, affecting computations of the 400d survey 
volume coverage; this uncertainty is mostly measurement in 
nature because we derive the Lx-M^oi relation internally from 
the same cluster set. Its effect is negligible for the high-M end 
of the mass function and becomes comparable to Poisson er- 
rors for low-M clusters. A representative compilation of the 
effects of -Mtot uncertain ties on the V{M) function is pre- 



sented in Fig. 15 of Paper II 

The general reliability ot our analysis is greatly enhanced by 
using independent, high-quality X-ray indicators of the total 
cluster mass — Tx, Mgas, ^x- Since the masses estimated from 
these proxies depend differently on the distance to the ob- 
ject, the high-z mass functions estimated with different prox- 
ies should agree only if the assumed background cosmology is 
correct. In principle, this can be used as an additional source of 
information for the distance-redshift relation and folded into 
the cosmological fit. However, this method is nearly equivalent 
to the /gas (z) test, which is more reUably carried which is more 
relibaly carried out using direct h ydrostatic mass estimates in 
relaxed clusters I Allen et al.|2008 1, and therefore we ignore this 
information. Instead, we use the agreement between different 
proxies observed for the best-fit cosmology as a comforting in- 
dication that there are no serious errors in our results. 

3. SUMMARY OF THEORY 

In the current paradigm of structure formation, galaxy clus- 
ters form via gravitational collapse of matter around large 
peaks i n the primordial density field (Kaiser 1984 Bardeen] 
et al.|198 6l. Their abundance and spatial distribution in a co- 



moving volume will thus depend on the statistical properties 
of the initial density field, such as gaussianity^ and power spec- 
trum (and hence the cosmological parameters that determine 
it), and could depend on the details of non-linear ampUtifi- 
cation of the density perturbations by gravity. Indeed, semi- 
analytic models based on the Unear primordial density field 
and a simple ansatz describing non-linear gravi tational col 



lapse of density peaks (Press & Schechter 1 974| Bond et al. 
1991 Lee & Shandarin| |1998[|Sheth et aL||20(yT) have proven 
to be quite successf ul m describing results of direct cosmolog- 
ical simulations of structure formation (e. 



, Lee & Shandarin 

,1999),Sheth et al. 200T] [Jenkins et al.,200T) . Ihe accuracy of 



the existing models, however, is limited and over the last sev- 
eral years the abundance of collapsed objects was calibrated by 
fitting appropriate fitting f unction to the results of direct cos- 
mologic al simulations ( Jenkins et al.'2001' 'Evrard et al."2002j 



Warren et al. 20061. The htting functions are expressed in the 
so-called universal form'" as a function of the variance of the 
density field on the mass scale M. The fact that such universal 
expressions exist implies that there is a direct link between the 
the Unearly evolving density field and cluster abundance. 

In our analysis we use the most recent accurate calibration of 
the halo mass function by Tink er et al.[(2008| , which provides 
fitting formulas for halo abundance as a function of mass, de- 
fined in spherical apertures enclosing overdensities similar to 
the mass we derive from ob servational proxies for the observed 
clusters. The [Tinker et al. fitting formulas are formally accu- 
rate to better than 5%Tor the cosmologies close to the concor- 
dance ACDM cosmology and for the mass and redshift range 

' We note however, that the current constraints on non-gaussianity from 
the CMB anisotropy measur ements imply that the expected effects on clusters 
are small jGrossi et al.|2 007|. 

In the sense tnat tne safne function and parameters could be used to pre- 
dict halo abundance for different redshifts and cosmologies. 
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Fig. 2. — Illustration of sensitivity of the cluster mass function to the cosmological model. In the left panel, we show the measured mass function and predicted 
models (with only the overall normalization at z = adjusted) computed for a cosmology which is close to our best-fit model. The low-z mass function is reproduced 
from Fig.[T] which for the high-z cluster we show only the most distant subsample (z > 0.55) to better illustrate the effects. In the right panel, both the data and the 
models are computed for a cosmology with Ha = 0. Both the model and the data at high redshifts are changed relative to the Oa = 0.75 case. The measured mass 
function is changed because it is derived for a different distance-redshift relation. The model is changed because the predicted growth of structure and overdensity 
thresholds corresponding to A^rit = 500 are different. When the overall model normalization is adjusted to the low-z mass function, the predicted number density 
of z > 0.55 clusters is in strong disagreement with the data, and therefore this combination of flm and Oa can be rejected. 



of interest in our study; at this level, the theoretical uncertain- 
ties in the mass function do not contribute significantly to the 
systematic error budget. Although the formula has been caU- 
brated using dissipationless N-body simulations (i.e. without 
effects of baryons), the expected effect of the internal redistri- 
bution of mass during baryon dissipation on halo mass func- 
tion are expected to be < 5% ( Rudd et al. 2008 ) for a realistic 
fraction of ba ryons that condense s to form galaxies . 

Similarly to | Jenkins et"ari l |2001| l and |Warren et al.| ( |2006[ >, the 
Tinker et al. formulas tor the halo mass function are presented 
as a function of variance of the density field on a mass scale M. 
The variance, in turn, depends on the Unear power spectrum of 
the cosmological model, P(fc), which we calculate as a product 
of the initial power law spectrum, k", and the transfer func- 
tion for the given mixture of CDM and baryons, computed 
using the analytic approximations of Eisenstein & Hu ( 1999| l. 
This analytic approximation is accurate to better than 2% tor 
a wide range of cosmologies, including cosmologies with non- 
negUgible neutrino contributions to the total matter density. 

Our default analysis assumes that neutrinos have a negUgi- 
bly small mass. The only component of our analysis that could 
be affected by this assumption is when we contrast the low- 
redshift value of derived from clusters with the CMB power 
spectrum normaUzation. This comparison uses evolution of 
purely CDM+baryons power spectra. The presence of Ught 
neutrinos affects the power spectrum at cluster scales; in terms 
of ffg, the effect is roughly proportional to the total neutrino 
density, and is 20% for Y,^v = 0.5 eV (we calculate the ef- 



4.2.8 of |Dunkley et al. 2008| l and so we experime nt al so with 
eutrino masses outside the Ly-a forest bounds (§|8.5ll. 



feet of neut rinos using the transfer function model of Eisen- 
stein & Hu 1999 1. Stringent upper limits on the neutrino mass 
were reported trom comparison of the WMAP and Ly-a forest 
data, X; Wv < 0. 17 eV at 95% CL ( Seljak et al. 2006 1. If neutrino 
masses are indeed this low, they would have no effect on our 
analysis. However, possible issues with modeling of the Ly-a 
data have been noted in the Uterature (see, e.g., discussion in 



4. FITTING PROCEDURE 



We obtain parameter constraints using the likeUhood func- 
tion computed on a full grid of cosmological parameters affect- 
ing cluster observables (and also those for external datasets). 
The relevant parameters for the cluster data are those that affect 
the distance-redshift relation, as well as the growth and power 
spectrum oflinear density perturbations: flM> ^a, w (dark en- 
ergy equation of state parameter), Cg (Unear ampUtude of den- 
sity perturbations at the 8 h^^ Mpc scale at z = 0), h, tilt of the 
primordial fluctuations power spectrum, and potentially, the 
non-zero rest mass of light neutrinos. This is computationally 
demanding and we describe our approach below. 

The computation of the Ukelihood function for a single com- 
bination of paramete rs is relat ively straightforward. Our pro- 
cedure (described in Paper ll\ uses the full information con- 
tained in the dataset, without any binning in mass or redshift, 
takes into account the scatter in the Mtot vs. proxy relations 
and measurement errors, and so on. We should note, how- 
ever, that since the measurement of the Mgas and Yx proxies 
depends on the assumed distance to the cluster, the mass func- 
tions must be re-derived for each new combination of the cos- 
mological parameters that affect the distance-redshift relation 
— D.My w, Oa, etc. Variations of h lead to trivial rescalings of 
the mass function and do not require re- computing the mass 
estimates. Computation of the survey volume uses a model for 
the evolving Lx - Mtot relation (see § 5 in Paper II), which is 
measured internally from the data and thus also depends on 
the assumed d{z) fiinction. Therefore, we refit the Lx - Mtot 
relation for each new cosmology and recompute V{M). Sen- 
sitivity of the derived mass function to the background cos- 
mology is illustrated in Fig.[2] The entire procedure, although 
equivalent to full reanalysis of the Chandra and ROSAT data. 
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TABLE 1 

COSMOLOGICAL CONSTRAINTS FROM X-RAY CLUSTER DATA 



Parameter Value Determined by Systematic errors Dominant source 

of systematic uncertainties 

D.Mh 0. 184 ± 0.024 Shape of the local mass function, §[5] ±0.027 Slope of the L - M relation. 

CIm 0.255 ± 0.043 Shape of the local mass function pfils ±0.037 Slope of the L~ M relation. 

HST prior on h, §[5] 

(T8(nM/0.25)°-^' 0.813 ±0.013 Normalization of the local mass ±0.024 Absolute mass calibration at z = 0. 

function, §[6| 

D.M 0.34 ±0.08 Evolution oTthe 73( -based mass ±0.055 Evolution of the M - T relation 

function, §[7] 



can be organized very efficiently if one stores the derived Pg{r) 
and r(r) computed in some reference cosmology. It takes 
Ri 20 sec on a single CPU to re-estimate all masses, refit the 
Lx - Mtot relation, and recompute volumes for each new com- 
bination of the cosmological parameters. 

The next step is to compute, for each combination of D.m, 
Oa etc., the likelihood function on a grid of those parameters 
which do not affect the distance-redshift relation. In our case, 
these are erg. h, and when required, the power spectrum tilt or 
neutrino mass. The cluster datad are extremely sensitive to erg 
and so we need a fine grid for this parameter. Fortunately, the 
mass function codes compute the mass functions for different 
values of as with other parameters fixed at almost no extra ex- 
pense. The sensitivity of the cluster data to h and tilt is much 
weaker, therefore the Ukelihood can be computed on a coarse 
grid for these parameters and then interpolated. 

With the acceleration strategies outlined above, it took us 
~ 9600 CPU-hours (or 20 days using multiple workstations) 
to compute the cluster Ukelihood functions on full parame- 
ter grids for several generic models (non-flat ACDM, constant 
dark energy equation of state in a flat universe, constant w with 
non-zero neutrino mass, linearly evolving w in flat universe, 
constant w in no n-flat universe). A lternatively, simulating the 
Markov chains (Lewis & Bridle 2002) with sufficient statistics 
for all these cases would require approximately the same com- 
puting time. 

After the cluster Ukelihood function was computed, we also 
computed for external cosmological datasets — WMAP 
(5-year results), Baryonic Acoustic Oscillations, and Super- 
novae la bolometric distances. Since we basically use ana- 
lytic Gaussian priors for these datasets (see § 8.1 below), these 
computations are fast and can be made on a hne parameter 
grid. We also use a Gaussian prior for the Hubble constant, 
h = 0.72 ±0.08, b ased on the results from the HST Key Project 
(jFreedman et al. 200 I I. This prior is important only when the 
constraints from the shape of the mass function (§[5| come into 
play and when external cosmological datasets are not used in 
the constraints. When fitting the cluster data, we also keep 
the absolute barvon density fixed at t he best-fit WMAP value, 
n^h^ = 0.0227 (Dun kley et al.|2008| l. This parameter slightly 
affects the c alculation of the Unear power spectrum (Eisenstein 
|&Hu|1998^ . When we add the WMAP information to the total 
constraints, we marginalize the WMAP likeUhood component 
over this parameter. If not stated otherwise, our cosmological 
fits also assume a pr imordial density fl uctuation power spec- 



trum with « = 0.95 (Spergel et al. 
pletely insensitive to variations o: 



.|2007 1. Our results are com- 
n within the WMAP mea- 
surement uncertainties and even to setting n = I. 

Once the combined likeUhood as a ftinction of cosmological 
parameters is available, we use the quantity -2 InL, who se sta 



In addition to statistical uncertainties, we also consider dif- 
ferent sources of systematics. We do not include systematic 
errors in the likeUhood ftinction but instead refit parameters 
with the relations affected by systematics varied within the es- 
timated Icr uncertainties. This approach allows as not only 
to estimate how the confidence intervals are expanded from 
combination of all systematic errors, but also to track the most 
important source of uncertainty for each case. A full analysis 
of systematic errors is presented in § |8.4| for the case of con- 
straints on constant w in a flat universe; in other cases the sys- 
tematic uncertainties contribute approximately the same frac- 
tion of the total error budget. We also verified that in the con- 
stant w case, our method of estimating the systematic errors 
produces the results which are very close to the more accurate 
procedure using the Markov chain analysis. 

5. CONSTRAINTS FROM THE SHAPE OF THE LOCAL 
MASS FUNCTION: Quh 

The shape of the cluster mass function reflects the shape of 
the Unear power spectrum in the relevant range of scales, ap- 
proximately 10 h Mpc in our case. This shape, for a reason- 
able range of parameters in the CDM cosmology is controlled 
( Bardeen et al. 1 9861 mostly by the quantity Hm^- It is useful 
to consider constraints on this combination separately because 
they are nearly independent of the rest of the cosmological pa- 
rameters we are trying to measure with the cluster data. 

Fixing the primordial power spectrum index to the WMAP 
value, n = 0.95, the fit to the local mass function" gives 
Cluh = 0.184±0.024 (purely statistical 68% CL uncertainties). 
The best fit value is degenerate with the assumed primordial 
power spectrum index, and the variation approximately fol- 
lows the relation AQm^ = -0.31A«. The variations of n within 
the range constrained by the WMAP data, ±0.015, lead to neg- 
Ugibly small changes in our derived Ciuh. 

An additional source of statistical uncertainty is that related 
to the derivation of the L - M relation, since we derive this 
relation from the same set of clusters. Uncertainties in the 
L - M relation are translated into those of the survey volume 
and hence the cluster mass function. Most of our cosmological 
constraints are primarily sensitive to the cluster number den- 
sity near the median mass of the sample. This median mass, 
the V(M) un certainties are small compared to statistics (see 
§6 in Paper II| l. The Cluh determination, however, is based on 



tistical properties are equivalent to the x distribution ( Cash 
1979 1, to find the best fit parameters and confidence intervals. 



' ' Including the high-redshift data, we obtain a consistent value, Cm 'i = 
0. 198±0.022. Combined with the HST prior on h, this leads to a measurement 
of Cm = 0.275 ± 0.043. However, using the high-z data makes the Cluh 
constraints dependent on the background cosmology and therefore we prefer 
to base this measurement only on the local mass function. Also, we use the 
Yx-based mass estimates for this and dg analyses. The other observables, Tx 
or Mgas, give essentially identical results, because a ll of them were normalized 
using the same set of low-z clusters (see jPaper II| for details). The difference 
between mass proxies is only important lor ttie measurements based on the 
evolution of the high-z mass function (§[7j. 
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Fig. 3. — Constraints on the dg and parameters in a flat ACDM cosmol- 
ogy from the total (both low and high-redshift) cluster sample. The inner solid 
region corresponds to -2 A In L = 1 from the best-fit model (indicates the 68% 
CL intervals for one interesting parameter, see footnotelTs} and the solid con- 
tour shows the one-parameter 95% CL region (-2AlnL = 4). The dashed 
contour shows how the inner solid confidence region is modified if the nor- 
malization of the absolute cluster mass vs. observable relations is changed by 
+9% (our estimate of the systematic errors). 



the relative number density of dusters near the high and low 
mass ends of the sample. Since the volume is a fast-decreasing 
ftinction at low Ms, the V{M) variations are important. The 
most important parameter of the L - M relation in our case is 
the power law slope, a (see eq. 20 in Paper III. Variations of a 
within the errorbars (±0.14) ofthe best ht value lead to changes 
in the derived Om h of ±0.027. Adding this in quadrature to the 
formal statistical errors quoted above, we obtain a total uncer- 
tainty of ±0.035. We have verified that other sources of sys- 
tematics in the determination are much less important 
than those related to the L - M relation. 

In principle, a non-zero mass of light neutrinos has some ef- 
fect on the perturbation power spectrum at low redshifts. We 
checked, however, that their effect on the shape of the cluster 
mass function is negligibl e for any Y,^v wit hin the range al- 
lowed by the CMB data ( Komatsu et al. 2008| l. Therefore, neu- 
trinos do not affect our results on LIm^- 

Our determination of Cl^h = 0.184 ± 0.035 compares well 
with the previous measurements using cluster data and galaxy 
power spectra. Of th e previous cluste r results especially note- 
worthy is the work of Schuecke r et al.|p00 3 1 whose constraints 
are based not only on the shape of the mass function but also 
on the clustering of low-z clusters. Their value is = 
0.239 ± 0.056 (errors dominated by uncertainties in the con- 
version of cluster X-ray luminosities into mass; this source of 
uncertainty is avoided in our work by using high-quality X- 
ray mass proxies). Ciuh is measured accurately also by galaxy 
redshift surveys. The results from the 2dF and SDSS su rveys 
are Quh = 0.178 ± 0.016 and .223 ± 0.023, respectively dColel 

we rescaled to n 



1) 



0.95 their 
The individual errorbars 



|et al. 2005 Teg mark et al [2004 
Dest fit values reported tor n 
in galaxy survey results are smaller than those from the clus- 
ter data; however, a recent work by Percival et al. ' (2007c I sug- 
gests that the previous galaxy redshitt results may be aftected 




Fig. 4. — Comparison with other ag measurements. Solid region is our 68% 
CL region reproduced from Fig.[3|(this and all other confidence regions corre- 
spond to A;f = 1, see footnote] fijon page!?). Blu e contours sho w the WMAP 
3 and 5-year results from Spergel et al. 1 2Cfn7 1 and Dunkley et al. 1 2008 1 (dot- 
ted and solid contours, respectiveiy). For other measurements, we stiow the 
general direction of degeneracy as a solid line and a 68% uncertainty in ug at a 
representative value of D^- Filled circles show the weak lensing shear results 
from Hoekstra et al. 1 2006 1 and Fu et al. 1 2008 1 (dashed and solid lines, respec- 
tively). Open circle shows results from a cluster sample with galaxy dynamics 
mass measurements I Rines et al. 2007 1. Finally, open square shows the results 
from Reiprich & Bohrmger 1 20U2 approximately the lower bound of recently 
published X-ray cluster measurements). 

by scale-dependent biases on large scales. Indeed, there is a 
tension between the SDSS and 2dF values at ~ 90% CL and the 
difference is comparable to the errorbars of our measurement. 

The cluster results can be improved in the future by extend- 
ing the range of the mass function measurements. Not only 
can this improve statistical errors in the mass function mea- 
surements but it can also improve the accuracy of the L - M re- 
lation, a significant source of uncertainty in our case. We note 
that it is more advantageous to increase statistics in the high-M 
range than to extend the mass function into the galaxy group 
regime. In addition to greater reliabihty of the X-ray mass esti- 
mates in the high-M systems, the surveys become dominated 
by cosmic variance approximately below the lower m ass cut in 



our sample (the cosmic var iance is estimated in §7.1 of Paper II 



using the prescription of Hu & Kravtsov|2003| l. 

Combined with the HSl' prior on the Hubble constant, our 
constraint on Ciuh becomes a measurement for the matter 
density parameter, Hm = 0.255 ± 0.043 (stat) ±0.037 (sys), 
where systematic errors are also dominated by the slope of 
the L - M relation. This agrees within the errors with other 
independent determinations, such as a combinat ion of BAO 
and CMB acoustic scales, Om = 0.256 ± 0.027 (P ercival et al.| 
|2007b ), and a combination of gas fraction measurements in 
massive clusters with the average baryon density from Big 
Bang Nucleosynthesis, Om = 0.28 ± 0.06 l |Allen et al.|2008i l. 
It also agrees with another independent nieasurement baseid 
on our data, Hm = 0.30 ± 0.05 from evolution of the cluster 
temperature function, see (§[7]below). 

6. CONSTRAINTS FROM THE NORMALIZATION OF 
THE CLUSTER MASS FUNCTION: Cg - D,m 
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The normalization of the cluster mass function is exponen- 
tially sensitive to a^, the amplitude of linear perturbations at 
the length scale Sh"^ Mpc, ap proximately corresponding to 
the cluster mass scale IjFrenk et al. 19901. Measuring this pa 
rameter with the cluster dSa has been a popular topi c of re 



_search, especially using statistics of X-ray clus ters (Prenk e t al. 

:etar 



|T990; Henry &Arnaud 199L Lilje 1992 Whitej 



|1993| and 

many others thereafter). The strong sensitivity of the predicted 
cluster number density to Og makes the determination of this 
parameter relatively insensitive to the details of the sample se- 
lection. Historically, different studies using very different clus- 
ter catalogs yielded similar results, if the data were analyzed 
uniformly. Determination of Cg is more s ensitive to calibration 
of the absolute mass scale. For example Pierpaoli et al. 1 2003 1 
show that if Mtot for a fixed value of Tx is varied by a factor 
of 1.5, ffg derived from the local cluster temperature function 
is changed by Acg r; 0.13. Smaller biases are introduced if 
the effects of scatter in deriving the mass-luminosity relation 
are neglecte d resulting in incor rect computations of the sur- 
vey volume I Stanek et al.|2006^ . Our present work includes 
advances in both of these areas and thus it is worth presenting 
an updated measurement of erg. 

Determination of Cg from the cluster abundance data usu- 
ally shows a strong de generacy with the Om parameter, typi 



cally ffg cx a J" (e.g., Huterer & White 2002 1. The nature of 



this degeneracy is that the mass function determines the rms 
amplitude of fluctuations at the given Mtot scale. The corre- 
sponding length scale is a function of Om (M ~ Om^^) and 
thus the derived Cg depends also on Hm and more weakly on 
th e local slope of th e linear power spectrum (see discussion 
in [White et aL|1993| We need, therefore, to constrain Cg and 
Om jointly. We used a grid of parameters of the flat'^ ACDM 
model (Hm, h, ffg), and computed the cluster likelihood us- 
ing the mass function for the local sample. We then add the 
Hubble constant prior (§[4}, and marginalized the combined 
likelihood over h. 

The results are shown in Figjsj'. For a fixed Clu, the value 
of ffg is constrained to within ±0.012 (statistical). The de- 
generacy between ffg and Om can be accurately described as 
ffg = 0.813(aM/0.25)"° *^. The Om range along this line is 
constrained by the shape of the local mass function combined 
with the HST prior on the Hubble constant (§|5}. Including the 
high-redshift data, we obtain very similar results. For example, 
forflM = 0.25, the total sample gives ffg = 0.803 ±0.0105, to be 
compared to ffg = 0.813 ± 0.012 from low-z clusters only. This 
implies that the ffg measurement is dominated by the more ac- 
curate local cluster data, as expected. 

Systematic errors of the ffg measurement are dominated by 
the uncertainties in the absolute mass calibration. To test the 
effect of these uncertainties, we changed the normalization of 
the mass vs. proxy relations by ±9% (our estimate of systematic 
errors in the mass scale calibration, see §|2]l. The effect, shown 

" The assumption of flatness (and background cosmology in general) has 
a minor effect on determination of ag because the measurement is dominated 
by the low-redshift sample. However, we note that when we use the ug infor- 
mation in the dark energy constraints (§[8]and thereafter), we do not use the 
results from this section directly. When we fit w, is effectively re-measured 
from the cluster data for each background cosmology. 

" The contours in these and subsequent figures correspond to the 95% CL 
region for one interesting parameter (A)(^ = 4). The inner solid region cor- 
responds to A)(^ = 1. This choice is made to facilitate quick estimates of the 
single-parameter uncertainty intervals directly from the plots. The total extent 
of the A^^ = 1 region in either direction is a good estimate for the 1-parameter 
68% CL interval (Cash,1976| L Similarly, the width of this region is a 68% CL 
interval assuming ttiat me second parameter is fixed. 



by the dotted contour in Fig.jsj is to shift the estimated values 
of ffg by ±0.02, just outside the statistical 68% CL uncertainties. 
This range can be considered as a systematic uncertainty in our 
ffg determination for a fixed Om- 

Our cluster constraints on ffg are more accurate (for a fixed 
Om) than any other method, even including systematic errors 
(Fig.|4j. It is encouraging that our results are in very good 
agreement with recent results from other methods. The mea- 
surements based on lensing sheer surveys, cluster mass func- 
tion with Mtot estimated from galaxy dynamics, and WMAP 
(5-year results assuming flat ACDM cosmology) are all within 
their respective 68% CL uncertainties from our best fit. This 
independently confirms that our calibration of the cluster mass 
scale is not strongly biased. Furthermore, the present system- 
atic errors in the cluster analysis are smaller than the statistical 
accuracy provided by WMAP-5 and other methods. This al- 
lows us to effectively use the ffg info rmation in the dark energy 
equation of state constraints (§ 8.3 1. 

We now move to models where the crucial role is played by 
the high-redshift cluster mass fiinction data. The first case to 
consider is combined constraints for Om and Oa in the non- 
flat ACDM cosmology. To better demonstrate what role the 
different components of the information provided by the clus- 
ter mass function play in the combined constraints, we con- 
sider two cases: a) when the full cluster mass function infor- 
mation is used, and b) when the shape information is artifi- 
cially removed thus leaving only the evolutionary information. 

7. CONSTRAINTS FOR NON-FLAT ACDM 
COSMOLOGY: Om - Oa 

In the first case, for each combination of parameters, we 
compute the full likelihood for the low and high-z mass func- 
tions and add the HST prior on the Hubble constant (this 
is necessary for effective use of the mass function shape in- 
formation, see §[4] and §[5|. We then marginalize the com- 
bined likelihood over non-essential parameters (ffg and h in 
this case) keeping the primordial power spectrum index fixed 
at the WMAP best-fit value, n = 0.95. Removal of the shape 
information (our second case) is achieved by letting n vary and 
marginalizing over it. This is approximately equivalent to us- 
ing a free shape parameter for the CDM power s pectra, the ap 
proch often used in earlier cluster studies (e.g., 
|2001| l. Constraints for both cases were obtained : 
tions estimated using all our three proxies, Tx 



Borgani et al. 
or mass tunc- 

■ ...-Mgas, and Yx- 

The results are presented in Figs![5]and[6] First, we can eas- 
ily identify the role of using the mass function shape informa- 
tion (illustrated for the Mgas and Yx proxies). Clearly, it mostly 
breaks the degeneracies along the Om axis. The best fit values 
and statistical uncertainties for Om are very close to those de- 
rived from the shape of the local mass function (and nearly 
identical to those from the total sample, §[5|. 

For a fixed Om, the observed evolution in the cluster mass 
function provides a constraint on Oa. Degeneracies in the 
Om - Oa plane provided by different mass proxies applied to 
the same set of clusters differ because of the different distance 
dependence of the Mtot estimates via Tx, Mgas, and Yx (see be- 
low). Even without the shape information, evolution in the Yx 
and Mgas -based mass functions requires Oa > at the 85% 
and 99.7% CL, respectively. Including the shape information, 
we obtain Om = 0.28 ± 0.04, Oa = 0.78 ± 0.25 (and Oa > is 
required at the 99% CL) from the Yx -based analysis. The evo- 
lution of the Mgas -based mass function gives Om = 0.27±0.04, 
Oa = 0.83 ± 0.15, and Oa > at 99.98% CL. The Tx-based 
mass function does not strongly constrain Oa but provides an 



VIKHLININ ET AL. 



0.6 



< 

C! 



OA 



0.2 




evolution 



evolution Mj 



0.2 



0.4 



0.6 



Qm 



Fig. 5. — Constraints for non-flat ACDM cosmology from evolution of 
the cluster mass function. The results using only the evolution information 
(change in the number density of clusters between z = and z « 0.55) are 
shown in blue and green from the Mgas and Tx-h&std total mass estimates. 
The degeneracies in these cases are different because these proxies result in 
very different distance-dependence of the estimated masses (see text for de- 
tails). The constraints from the Yx-based mass function are between those 



for Mgas and Tx (Fig.|6j. Adding the shape of the mass function information 
breaks degeneracies with Om, significantly improving constraints from Mgas 
and Yx with little effect on the Tx results. 



independent measurement of Om with almost no degeneracy 
with Oa: Om = 0.34 ± 0.08, in good agreement with the mass 
function shape results (and also previous measurement s based 
on ev olution of the cluster temperature function, see Henry 
|2004| ). In a flat ACDM model (the one with Hm + = Ijthe 
constraint is sUghtly tighter, Clu = 0.30 ± 0.05. 

Systematic uncertainties of the D.a measurements are domi- 
nated by possible departures of evolution in the Mtot vs. proxy 
relations. This issue is discussed in detail below in connec- 
tion with the dark energy equation of state constraint (§ |8.4[ l; 
here we note only that the systematic uncertainties are approx- 
imately 50% of the purely statistical errorbars on the dark en- 
ergy parameters {D.a, w). Therefore, our cluster data provide 
a clear independent confirmation for non-zero Ha- 

Comments on the role of geometric information in the cluster mass 
function test— Cosmological constraints based on fitting the 
cluster mass function generally use not only information from 
growth of structure but also that from the distance-redshift re- 
lation because derivation of the high-z mass functions from 
the data assumes the d{z) and E{z) functions. Quite gener- 
ally, the estimated mass is a power law function of these de- 
pendencies, M oc d{zY E{zY^. Different mass proxies have 
different j3 and £, and thus combine the geometric and growth 
of structure information in different ways and lead to differ- 
ent degeneracies in the derived cosmological parameters. We 
find that strongly distance-dependent proxies (such as Mg^s, 
see Paper II I are intrinsically more powerful in constraining 
the dark energy parameters (Oa, w). By contrast, distance- 
independent proxies such as Tx result in poor sensitivity to 
dark energy but instead better constrain Om- This is well illus- 
trated by the results in Fig.js] The Mgas based estimates for Mjot 
result (if we ignore the shape of the mass function information) 
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Fig. 6. — Same as Fig.|5]but for Y^-based mass estimates. 



in degeneracy approximately along the line Om + ^^a = 1- In 
fact, the evolution of the cluster mass functions derived from 
Mgas can be made broadly consistent with the Hm ~ 1. ^^a ~ 
cosmology if one allows for strong deviations from the CDM- 
type initial power spectra ( Nuza & Blanchardl|200 6 1 . How- 
ever, the mass functions estimated from the temperatures of 
the same clusters are grossly inconsistent with such a cosmol- 
ogy, irrespective of the assumptions on the initial power spec- 
trum (Hm = 1 is 8.3(Tawayfrom the best fit to the temperature- 
based mass function, Fig.js}. It is encouraging that the 68% 
CL regions for all three mass proxies overlap near the "concor- 
dance" point at Om = 0.25 - 0.3 and Oa = 0.7 - 0.75. 

8. FLAT UNIVERSE WITH CONSTANT DARK ENERGY 
EQUATION OF STATE: Wq - Ox 

Next, we study constraints on a constant dark energy equa- 
tion of state, tvo = px/px, in a spatially flat universe. The anal- 
ysis using cluster data only is equivalent to the Om - ^^a case 
(§[7|. We compute the Ukelihood for the cluster mass func- 
tions on a grid of parameters: present dark energy density D.x 
(= 1 - flivi). Wo, h, and ffs, then add the HST prior on the Hub- 
ble constant (§[4|. MarginaUzation over non-essential parame- 
ters, h and Cg. gives the Ukelihood as a function of Hm and Wq. 
We also obtain the equation of state constraints combining our 
cluster data with the thr ee external cosmolog ical data sets (fol- 

for the choice of 



lowing the reasoning of Dunkley et al.|2008 

these datasets) 



8.1. External Cosmological Datasets 

SNIa — We use the distance moduU estimated for the Type la 
supern ovae from t he HST sample of Riess et al. (20071, SNLS 
survey ( Astier et al. 2006 1, and ESSENCE survey (Wood-Vasey| 
et al.|2007 l, combined with the nearby supernova sample (we 
used a combination of all these samples compiled by Davis] 
et al.|2007) l. Calculation of the SN la component of the like- 
iihood function for the given cosmological model is standard 
and can be found in any of the above references. 
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TABLE 2 

Parameter constraints from combination of clusters with other cosmological datasets 





Parameter Value 


Dataset 


Systematic 
errors 


Dominant source 
of systematics 


Wo . 


Flat (Oj. = OJ, constant w (w = woJ 
-1.14±0.21 


evol+shape+ft, §|8.2| 


±0.10, ±0.08, 


Evolution of Mtot vs. proxy relations, 
evolution in Lx - Mtot, respectively 


Wo . 
Wo . 


-1.08 ±0.15 

-0.97 ±0.12 


evol+cmb, §|8.3 
evoI+cmb+(78+ 


5ao 




±0.025 
±0.038 


Evolution of Mtot vs. proxy relations 
Absolute cluster mass calibration 


Ox 
Wo . 


0.732 ±0.016 

-0.995 ± 0.067 


cmb+bao+sn, § 
cmb+bao+sn, § 


8.3 




±0.076 


SN systematics 



Wo . 

fix 
h.. 

0-8 . . 



-0.991 ± 0.045 
0.740 ± 0.012 
0.715 ±0.012 
0.786 ± 0.011 



cmb+bao+sn+evol+(J8, § 
cmb+bao+sn+evol+ijg, § 
cmb+bao+sn+evol+(T8, § 
cmb+bao+sn+evol+(T8, § 



±0.022, ±0.033 SN systematics, cluster masses 



Flat (nj. = OJ, constant w (w = wo), non-zero neutrino mass 

Wo -1.02 ±0.055 (T8+cmbv+cmb+bao+sn+evol, § 

Y.my 0.1 ± 0.12 eV, (T8+cmbv+cmb+bao+sn+evol, § 

< 0.33 eV (95% CL) 



■K3 



±0.064 SN systematics. 

±0.1 eV SN systematics, cluster masses 



Flat ("nj- = OJ, evolving w: w = wo + iVo(l-a) 

Wa + 3.64(1 + wq) 0.05 ±0.17 cmb+bao+sn+evol+(T8, §|9.1| 

Constant w (w = wo), non-flat (Cli^ + 0) 
-1.03 ±0.06 



Wo 



cmb+bao+sn+evol+(J8, § 9.2 



Note. — Codes used in column 3: evol — evolution of the cluster mass function; h — HST prior on Hubble constant; shape — shape of the cluster mass 
function; cmb — WMAP-5 distance priors; og — comparison of the cluster-derived og with the CMB power spectrum normalization (reflecting growth of 
perturb atio ns between zcmb and z = 0); bao — BAO distance prior; sn — SN la luminosity distances; cmby — WMAP-5+BAO+SN constraints on neutrino 
mass (§[8!5). 



Baryonic Acoustic Oscillations — Detection of the barynic acous- 
tic peak in the correlation function for large red galaxies in the 
SDSS survey leads to a good measurement of the combination 



{czYH{z) 



-,1/3 



0.98 



0.469 ±0.017 (1) 



at z = 0.35 ( Eisenstein et al.|2005 "SDSS LRG sample"). This 
prior mostly constrains Om but has some sensitivity also to the 
dark energy equation of state. 
A more recent measurement of the BAO pea ks in the com- 



bined SDSS and 2dF survey data is presented in Percival et al. 
(|2007ajl who determine the BAO distance measiire at two red^ 
shifts (z = 0.2 and z = 0.35) instead of one in [ Eisenstein 
2005 1. These new data are somewhat in tension { 

Fig 



let al. 



2a) 



with the SN-i-WMAP results (see, e.g., Fig. 1 1 in Percival et al. 
l2007al, which may artificially tighten the constraints when the 
BAO data are combined with SN la, WMAP, and clusters. We 
checked, however, that from the combination of SN la, WMAP, 
and SDSS-LRG BAO, we derive the paramete r constraints tha t 
are essentiall y equivalent to those in Komatsu et al.| (|2008|, 
who used the [Percival et al. priors. Therefore, the choice of 
the BAO dataset is unimportant in the combined constraints. 

WMAP-5— The likelihood for WMAP 5-year data is c om- 
puted using a si mphfied approach described in § 5.4 of Ko- 
matsu et"dr|(|200 8 l. This involves a computation, for a given 
set of cosmological parameters, of three CMB parameters — 
angular scale of the first acoustic peak, the so called shift 
parameter, R; and the recombination redshift, z». The like- 
lihood for the WMAP-5 data is then computed using the 
coyariance matrix for R, and z* provided in Komatsu 



et al. This method is almost as accurate as direct comp uta 
tion of the WMAP Ukelihood (Wang & Mukherjee'2007'i but 
is much faster, which allowed us to explore the entire multi- 



dimensional grid of the cosmological parameters instead of 
running Markov chain simulations. One additional note is that 
to compute the CMB likelihood, we had to add the absolute 
baryon density, flfoft^, to our usual set of cosmological param- 
eters and then marginalize over it. The reason is that while 
the average baryon density has very little impact on the rest 
of our analysis, the CMB data are very sensitive to 0.i,h^, thus 
any variation of h must be accompanied by the corresponding 
variation of fl;, without which the computation of the CMB 
Ukelihood would be inadequate. 

The method outUned above recovers essentially the entire 
information from the location and relative amplitudes of th e 
peaks in the CMB power spectrum ( Wang & Mukherjee|2007| l. 
One additional piece of information is the absolute normal- 
ization of the CMB power spectra, reflecting the amplitude of 
density perturbations at the recombination redshift, z,t « 1090. 
Contrasted with iTg determined from our cluster data at z « 0, 
it constrains the total growth of density perturbations between 
the CMB epoch and the present, and thus is a powerful addi- 
tional dark energy constraint. 



WMAP-5 plus local as — The WMAP team provides the am- 
plitude of the curvature perturbations at the k = 0.02 Mpc"^ 
scale. 



(2.21 ±0.09) X 10" 



(2) 



Section 5.5 in Komatsu et al. (2008| gives the prescription of 
how to predict this observable tor a given set of cosmological 
parameters and ffg- A useful accurate fitting formula can also 
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Fig. 7. — Constraints on the present dark energy density Clx and constant 
equation of state parameter wo derived from cluster mass function evolution 
in a spatially flat Universe. The results for Mgas and Ix^based total mass esti- 
mates are shown in red and blue, respectively. The inner solid red region shows 
the effect of adding the mass function shape information (|5j to the evolution 



be found in |Hu & Jain| ( |2004| l: 



1.79 X 10-* \ 0.024 / \ 0.14 / 

(, .-0.693 
072) 



(3) 



0.76 



(we adjusted n umerical coefficients to take into account that 
the Hu & Jain approximation uses the CMB amphtude at fc = 
0.05 Mpc"^ while the WMAP-5 results are reported for k = 
0.02 Mpc~^). In this equation. Go is the perturbation growth 
factor between the CMB redshift and the present, normal- 
ized to the growth function in the matter- dominated universe: 
G(z) = (1 + z) S{z)/S{zcmr)- This fitting formula helps to 
understand the nature of the ffs vs. CMB amphtude constraint. 
The relation between Cg and AfH depends on the absolute mat- 
ter and baryon densities, Om^^ and D.i,h^ (well -measured by 
the CMB data alone), and on the total growth factor. Go, and 
the absolute value of the Hubble constant, h. Both of these 
quantities provide powerful constraints on any parametriza- 
tion of the dark energy equation of state (Hu 200 51, and their 
combination does so as well. 

Inclusion of this information in the total likelihood is 
straightforward. Given the usual set of cosmological param- 
eters {D.X, Wo. h) plus (Ts. one computes 



(A^ X 10' -2.21)^0.09^ 



(4) 



XcMBnorm 

where Atn c an be obtained either from eq. js] or as described in 
Komatsu et al. (2008 1. The XcMBnorm component is then added 
to the cluster and the sum marginalized over erg. 

8.2. Wofrom Cluster Data Only 

Constraints on the present dark energy density Clx and con- 
stant equation of state are presented in Fig.|7] For compar- 




0.85 



Fig. 8. — Comparison of the dark energy constraints from X-ray clusters and 
from other individual methods (supernovae, baryonic acoustic oscillations, 
and WMAP). 

ison, we show separately the results derived only from evo- 
lution of the Mgas and Yx -based mass functions, and the ef- 
fect of including the mass function shape information (§[7]de- 
scribes the procedure for removing shape information from 
the cluster Ukelihood function). We do not consider here the 
Tx based mass estimates because they provide little sensitivity 
to the dark energy parameters {§p)- Just like in the Hm - ^^a 
case, evolution of the Mgas and Yx-'based mass ftinctions con- 
strains different combinations of Wq and fix- The width of the 
confidence regions across the degeneracy direction is similar 
but the gas-based results are less incUned giving a httle more 
sensitivityto Wo for a fixed dark energy density — Awg = ±0.17 
from the Mgas-based functions and Awq = ±0.26 from Yx- 

Adding the mass function information combined with the 
HST prior on h breaks the degeneracy along the fix direction. 
For example, the elhpse in FigJTlshows the 68% CL region from 
fitting both the evolution anoshape of the Mgas-based mass 
ftinction. The one-parameter confidence intervals in this case 
are fix = 0.75 ±0.04 and ivo = -1.14±0.21. These results com- 
pare favorably with those from other individual methods — 
supernovae, BAO, WMAP (Fig.jsf, although the supernovae 
and CMB data provide tighter constraints on Wq for a fixed 
fix- The real strength of the cluster data is, however, when they 
are combined with the CMB and other cosmological datasets. 
The combined constraints are very similar for the Mgas and Yx - 
based cluster mass functions, and therefore we discuss only the 
former hereafter. 

8.3. Wo from the Combination of Clusters with Other Data 

First, we consider a combination of the cluster data with 
the WMAP distance priors (see § 5.4 in Komatsu et al.|2008^ . 
Cluster data bring information on growth of density pertur- 
bations and normaUzed distances in the z ^ 0.0 - 0.9 inter- 
val, and — weakly — on the parameter. Adding this in- 
formation reduces the WMAP-only uncertainties on Wo and 
fix approximately by a factor of 2 (dark blue region in Fig.|9}: 
Wo = -1.08 ± 0.15, Ox = 0.76 ± 0.04. 

A much more significant improvement of the constraints 
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Fig. 9. — Dark energy constraints in a flat universe from the combination of 
the CMB and cluster data (dark blue region). Adding the vs. CMB normal- 
ization information significantly improves constraints on wo for a fixed Clx 
(inner red region). 



arises from the ffg determination from low-redshift clusters 
(dark red region in Fig.[9l. Comparison of the local determi- 
nation of (Tg with the CMB normaUzation mostly provides a 
measurement of the total perturbation growth factor between 
ZcMB and the present. This depends more sensitively on ivq 
than the evolution of the cluster mass function because of first, 
larger redshift leverage, and second, because the perturbation 
amplitude at high z is measured more accurately by CMB than 
by 37 clusters from the 400d survey. 

Is it appropriate to use the vs. CMB normalization infor- 
mation in the dark energy constraints or does it require un- 
reasonable interpolation of the dark energy parametrization 
to high redshifts? We note in this regard that for any com- 
bination of the cosmological parameters in the vicinity of the 
"concordance" model, Wq ~ -1, fix = 0.25 - 0.3, the Uni- 
verse becomes matter- dominated and enters the deceleration 
stage by z - 1.5-2; the growth of perturbations is basically 
fixed after that at G(z) = 1. In other words, the CMB data 
can be used to safely predict the ampUtude of density pertur- 
bations at z = 1.5-2 almost independently of the exact dark 
energy properties. As long as it is appropriate to use a par- 
ticular dark energy parametrization in the z = - 2 interval, 
it is therefore appropriate to use the same model for the joint 
clusters+WMAP fit. 

By itself adding the ffg information does not significantly 
improve the Wo and D.x constraints (the total extent of the la 
confidence regions is similar to the WMAP-i-evolution case), 
but the confidence region becomes much more degenerate 
with fix (see inner red region in Fig.[9}, which increases the 
potential for improvement when we combine these results with 
other cosmological datasets, BAO and supernovae. 

The combined constraints from all four cosmological 
datasets are shown in Fig. 10 (inner dark red region). The 68% 
one-parameter confidence mtervals are fix = 0.740 ± 0.012 
and tvo = -0.991 ± 0.045. The importance of adding informa- 
tion from our cluster samples is illustrated by a factor of ~ 1.5 
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Fig. 10. — Dark energy constraints in flat universe from combination of all 
cosmological datasets. We find wo = -0.9 91 ± 0.045 (±0.04 systematic) and 
Qx = 0.740 ± 0.012, see Table|2]and § [83] 

WMAP+SN-t-BAO data alone: we obtain Wq = -0.995 ± 0.067 
without clusters (dark blue region in Fig.fTo] these results ar e 
essentially identical to those reported in Komatsu et al. 2008| l. 
Perhaps more importantly, including the cluster data also re- 
duces systematic uncertainties by a similar amount ( §|8.4| l. 

The best-fit values of the Hubble constant and Cg from the 
combination of all datasets are h = 0.715 ± 0.012 and ffg = 
0.786 ± 0.011. These values are within 68% confidence inter- 
vals of their determination by direct measurements (HST Key 
Project results for h and fitting the low-z cluster mass fiinction 
for ffg). The best-fit combination of the dark energy parame- 
ters is also within the la confidence regions for each individ- 
ual dataset included in the constraints (Fig. 101. Therefore, the 
best-fit cosmological model is di good fit to the data. In particu- 
lar. Fig. 17 from Paper II shows that the mass function models 
computed in the ACDM cosmology (tVo = -1) provide a very 
good description of the data. 

8.4. Systematic Uncertainties in the Wq Measurements 

We estimate the effect of known sources of systematics on 
the cosmological constraints by varying the corresponding in- 
dividual sets of data or internal relations (e.g., evolution in 
Lx-Mtot entering the survey volume computations) within the 
estimated Icr interval. We assume, optimistically, that the cur- 
rent WMAP and BAO data are free from significant system- 
atics (i.e., that they are smaller than statistical uncertainties), 
and consider systematic errors only in the SN la and cluster 
datasets. In most cases, a single source clearly dominates the 
systematic error budget for a particular measurement, so we 
report on only those dominant sources. 

The largest known source of systematic error in the SN la 
analysis is the correction for extinction in host galaxies and 
uncertainties in intrinsic colors of SN la (e.g., Fri eman et eJ!) 
[2008 1. As a measure of systematic uncertainty in the co mbined 

, quoted bylWood 



SN sample we use ±0.13 in tvp for fixed fix 
[Vasey et aT ( |2007| l. We implement these errors by computmg 
the SN likelihood in our experiments for (flx> + 0.13) and 



reduction of the measurement uncertainties with respect to the 



■ experiments 
(flx> Wo - 0.13) instead of (flx> Wq). 
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8.4.1. Main Sources of Cluster Sustematics 

The largest sources of systematic errors in the cluster anal- 
ysis are those in the normalization of the Mtot vs. proxy re- 
lations. They can be separated into two almost independent 
components: 1) how accurately is the absolute cluster mass 
scale established by X-ray hydrostatic M(ot estimates in the 
low-redshift clusters, and 2) how accurately can we predict 
evolution in the Mtot vs. proxy relations, i.e., the relative mass 
scale between low and high redshift clusters. The first compo- 
nent mainly affects the measurements and associated dark 
energy constraints, while the second component affects the 
results derived from using only evolution in the cluster mass 
function (those in Fig.[7j. Our estimates of the Mtot systemat- 
ics are discussed extensively in |Paper II For the absolute mass 
scale (Mtot for fixed Yx, Tx, or Mgas) at z r* 0, we estimate 
AMsys/M < 9% mainly from comparison of the X-ray and 
weak lensing mass estimates in representative samples. This 
source of error is implemented by changing the normalization 
of the Mtot vs. Yx, Mg^s, or Tx relations at z = by ±9%. For 
uncertainties in the evolution of the Mtot vs. proxy relations, 
we estimate AM/M « 5% at z = 0.5, mainly from comparison 
of the prediction of different models describing observed small 
deviations of the cluster scaling relations from self-similar pre- 
dictions, and from the magnitude of these deviations and cor- 
responding corrections we apply to the data. These uncertain- 
ties are implemented by multiplying the standard scaling rela- 
tions by factors of (1 + z)^"-'^. 

Comparable to the evolution in the Mtot vs. proxy relation 
are measurement uncertainties in the evolution factor for the 
Lx - Mtot relation. We do not use Lx to estimate the clus- 
ter masses, but the relation is required to compute the survey 
volume for the high-z sample. The resulting volume uncer- 
tainty depends on the mass scale, and can become comparable 
to the Poisson error for the comoving cluster number density 
(see § 5. 1.3 i n|Paper II( I. We tested how this influences the cos- 
mological fit by varying the parameters of the Lx - Mtot rela- 
tion within their measurement errors around the best fit [the 
evolution of Lx for fixed Mtot in our model is parametrized as 
E(zy and y is measured to ±0.33, see § 5.1.3 in Paper II| . 

Other sources of systematics in the cluster analysis (sum- 
marized in Paper III are negligible compared to those out- 
lined above. We verified also that uncertainties in the intrin- 
sic scatter in the Mtot -proxy relations are not important. The 
main reason is that in the dark energy constraints, we use 
high-quality mass proxies (Yx and Mgas), which should pro- 
vide mass estimates with small, 7-10% scatter. Variations of 
this scatter by up to ±50% with respect to the nominal values 
do not significantly change the best fit cosmological parame 



ters. T his conclusion is seemingly different from Lima & Hu 
(2005 I because in that paper, they consider proxies with larger 



scatter (the effect on the cosmological parameter constraints 
is proportional to scatter squared), and also they assumed that 
the normalizations in the Mtot vs. proxy relation are obtained 
from self-calibration while we use direct mass measurements 
for a well-observed subsample. 

The variations of the best-fit parameters due to the systemat- 
ics discussed above are reported in Table[2]along with the dom- 
inant source of error for each combination of cosmological 
datasets. For example, variations in the evolution of the Mtot - 
Mgas and Mtot - Yx relations affect the best fit to the cluster data 
only by Awq = ±0.1 , wh ile statistical uncertainties are ±0.2 to 
±0.3 for fixed D.x (§ |8.2| ; unless the systematics in this case are 
a factor of two larger than our estimates, they are unimportant. 



8.4.2. Systematics in the Combined Constraints 

The most interesting case to consider is reduction in the 
systematic errors from combining both SN and cluster data 
with the WMAP and BAO priors. In the SN+CMB+BAO case, 
the supernovae systematics cause variations in the best-fit Wq 
by ±0.076 (reduced from ±0.13 for the SN-only case mainly 
by including WMAP priors). Cluster systematics affects the 
Wo constraints from the clusters-i-WMAP-i-BAO combination 
by ±0.04 (dominated by the ±9% uncertainties in the abso- 
lute mass scale). The influence of both sources of error is sig- 
nificantly reduced in the combined constraints. We find that 
the best fit Wq from SN+clusters+WMAP+BAO is affected by 
±0.022 by SN systematics, and by ±0.033 by cluster systemat- 
ics. The total systematic error in the combined constraint is 
thus Awo = ±0.04, almost a factor of 2 reduction from ±0.076 
achievable without clusters. 

We also note that if we significantly underestimate the clus- 
ter systematics, the most likely direction is that the cluster total 
masses are underestimated". If cluster Mtot are revised high, 
this would lead to an increase in the derived erg. and decrease 
in Wo when cluster data are combined with the CMB priors. 
Dark energy models predicting the equation of state parame- 
ter significantly above Wo = - 1 will be even less consistent with 
observations in this case. 

8.4.3. Prospects for Futher Reduction of Systematic Errors 

It is reassuring that all sources of systematic errors we con- 
sidered affect the dark energy equation of state constraints 
within the statistical measurement errors. This implies that 
while systematic errors are important, they do not yet dom- 
inate the current error budget. The situation will reverse in 
the future as the datasets expand. More effort will be needed 
then to reduce the systematics still further We briefly out- 
line the prospects for reducing the cluster-related systemat- 
ics. Some of this will happen automatically as the high-z sur- 
veys become deeper and cover a larger area. For example, the 
V{M) uncertainties for our range of redshifts can be elimi- 
nated simply by decreasing the flux threshold by a factor of 
~ 4 compared to the 400d limit, making the sample volume- 
limited; such an extension will provide also a more accurate 
measurement of the Lx - Mtot relation. The absolute calibra- 
tion of Mtot in low-z clusters can be improved by constraining 
sources of non-thermal pressure (e.g., if turbulence is of any 
importance for the Mtot estimates, it is easily detectable with 
an X-ray microcalorimeter), or through stacked weak lens- 
ing analysis (e.g., measuring average lensing shear profiles for 
a large set of clusters with the same Yx). To improve limits 
on non-standard evolution in the Mtot vs. proxy relations, we 
cannot use direct mass measurements of the high-z objects 
because they will be degenerate with the assumed distance- 
redshift relation. Instead, we should improve reliability of nu- 
merical models for cluster evolution. The biggest uncertain- 
ties in these models at present are related to the processes of 
gas cooling and star formation, and also to energy feedback 
from the central AGN. The strategy for future progress can 
be based on the fact that these processes most strongly affect 
cluster cores, which we do not use for the mass estimates. We 

'* The X-ray hydrostatic analysis includes only the gas thermal pressure and 
assumes that the cluster gas body is close to being spherically symmetric. The 
presence of additional components in the pressure, dumpiness and turbulent 
motions in the gas all lead to underestimation of Mtot derived from X-ray data. 
Probably the only possibility for overestimation of Mtot in the X-ray analysis 
is a gross miscalibration of the Chandra spectral response, for which strong 
experimental limits are available. 
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Fig. 11. — Equation of states from WMAP, BAO, SN la, and clusters in the 
case of non-zero neutrino mass. 



can, therefore, use the data from the central regions to bracket 
a likely range of uncertainty in the model predictions for the 
cluster outer regions, where we derive the Mtot proxies. How- 
ever, even with the current estimated uncertainties, the sam- 
ples can grow by a factor of ~ 4 before the systematics start to 
dominate. Ultimately, as the cluster surveys detect ~ 10* clus- 
ters with accurately measur ed X-ray paramet ers, the so- called 
self-calibration techniques l |Majumdar & Moh r 2004; Lima & 



Hu|2004^ can be employed to further constrain the evolution 
mtnelWtot vs. proxy relations. 



itot vs. proxy 

8.5. Effects of Non-Zero Neutrino Mass 

If light neutrinos have masses in the range of a few 0.1 eV, 
they become non-relativistic between Zcmb and z = 0, and 
this transition produces distortions in the matter pertur- 
bations power spectrum relative to prediction of the pure 
CDM-^baryon s model. Using approxi mations of the transfer 
function from Eisenstein & Hu l |1999| l, it is easy to verify that 
the effect is approximately proportional to the total mass of 
neutrinos (more exactly, to Y, wJv/^^m), and the rms fluctua- 
tions at cluster scales today are suppressed by approximately 
20% if X fw,, = 0.5 eV and D.m = 0.26. This effect is far outside 
the measurement uncertainties in from clusters (we quote 
systematic errors of 3% from uncertainties in the Mtot cali- 
bration and statistical uncertainties are even smaller, see Ta- 
bleflj. Therefore, neutrino masses in this range a) may affect 
the aark energy constraints when cluster data are combined 
with WMAP (because they will effectively change the relation 
between Cg and the CMB normaUzation, eq.|3}, and b) can be 
tightly constrained by our cluster data. 

To test the effect of neutrinos, we ran an additional set of 
models in which the total neutrino mass was allowed to vary 
between and 1 eV. For simplicity we assumed that there are 
3 neutrino species with the same mass, but the final results are 
not very sensitive to this assumption. The only component 
of our procedure which is significantly affected by non-zero 
neutrino mass is contrasting the cluster- derived Cg with the 
WMAP normalization of the CMB power spectrum. We can 
no longer rely on eq.l[3]l and should instead use the full proce- 



dure described in § 5.5 of jKomatsu et al.|p008 1. Otherwise, the 
analysis is equivalent to the Y,mv = case. The Ukelihood for 
all cosmological datasets was computed on our usual grid plus 
2 as an additional free parameter, and then marginalized 
over D.xyh, and Cg . Finally, we took into account that a combi- 
nation of WMAP, BAO, and SN data provides some sensitivity 
to neutrino mass through the so-called e arly integr ated Sachs- 
Wolfe effect (see discussion in § 6.1.3 of Kom atsu et al.,2008l 



and references therein). From this analysis, Komatsu et al. de- 
rive a 95% upper limit of X n^ v < 0.66 eV. Since our procedure 
of using WMAP priors (§ |8.1[ l ignores this additional informa- 
tion, we included it approximately by adding a Gaussian prior 
= ± 0.33 eV to the final marginalized likelihood. 
The derived constraints on Y, and Wo are shown in 
Fig. 1 1 As expected, when the Cg vs. CMB normalization con- 
straint is added, there is a degeneracy between the best-fit Wq 
and total neutrino mass. If we were using only clusters and 
WMAP, the degeneracy would approximately follow the line 
Wo + 1 = -0.4 Y, 'Mv and would extend t oY^i^v ~ 1-3 eV (th e 
WMAP-only bound on neutrino mass, Dunkley et al.|200^ . 
This degeneracy is broken, however, when we add the B AO and 
SN information: low values of Wq required by clusters-i-CMB 
for high values of neutrino mass are inconsistent with these 
two datasets. Therefore, a combination of all four datasets can 
be used to constrain both Wg and neutrino mass. The best fit 
value is = 0.10 ± 0.12 eV, with a 95% CL upper limit 

of 2 < 0.33 eV. This limit is significantly tighter than that 
achievable without clusters (< 0.66 eV at 95% CL). The con- 
straint on Wo degrades somewhat compared to the = case: 
Wo = -1.02±0.055 (compared to ±0.045 for m^, = 0),butis still 
better than ±0.067 without clusters (see Table|2]l. To conclude, 
adding the cluster information allows us to set tight Umits on 
the neutrino mass while still improving the Wo measurements 
with respect to the SN-i-WMAP-i-BAO case. 

Our constraints on neutrino mass are still weaker than the 



published re sults from Ly-a forest data, Y'^v < 0.17 eV iSel- 
jak et aL(200 6). Both the cluster and Ly-a based constraints 
use the same effect — suppression of the power spectrum at 
small scales by neutrinos, — but they have completely differ- 
ent systematics. The main unknown in the Ly-a analysis is 
the thermal state of the low-density IGM, usually estimated 
from numerical simulations; it has been suggested that the 
thermal state may be more complex than assumed in previous 
work thus significantly weakening the m,, bounds (Bolton et al. 
2008| . For clusters, the main uncertainty is the absolute mass 
calibration for low-redshift objects which affects the measure- 
ment of (Tg (§|6j. The 9% systematic uncertainties on AM/M 
that we quote would translate into approximately ±0.075 eV 
for Y ffiy, negligible compared to the current statistical uncer- 
tainties. We note that if the X-ray cluster mass measurements 
are wrong by more than 9%, it is almost certainly in the sense 
that they are underestimated (see footnote 14 on page[T2}; the 



true value of Cg will then be higher than our measurement and 
the bound on neutrino mass will be even tighter. Therefore, 
our 95% CL bound oiY^v < 0.33 eV can be considered as a 
conservative upper limit. 

9. MORE GENERAL DARK ENERGY MODELS 

Finally, we demonstrate how our cluster data improves pa- 
rameter constraints for more general dark energy models. We 
consider two cases — evolving equation of state, w = w{z), and 
constant equation of state in a non-flat universe. The results are 
presented less completely than for the case of constant w in a 
flat universe. We also do not discuss systematic uncertainties 
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Fig. 12. — Constrains on evolving equation of state, w(z) = wo + Waz/{l+z), 
in flat universe. 

separately for these cases; we checked that the importance of 
different sources of systematics and their fraction of st atisti cal 
uncertainties is approximately the same as reported in § i 
the constant w, flat universe case. 



8.4 



for 



9.1. w{z) in Flat Universe 

We consider an often used parametrization of the equa- 
tion of state evolution in which w changes linearly with the 
expansion factor, w(fl) = Wo + Wfl(l-fl), or equivalently, 
w{z) = Wo + WflZ/(l + z). We do not consider more complex 
parametrizations because constraints on the evolution term 
are still weak, and because neither parametrization has a clear 
physical motivation. 

The likelihood function is computed on the Hm. Wq, Wa, h, 
(Ts grid and then marginalized over D.m, h, and Ug, leading 
to constraints in the Wq 
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Con- 



Wa plane shown in Fig. 
straints on Wa are weak with or without clusters. Fofexample, 
the model with Wq = -1.2 and Wa = I (leading to iv = -0.7 by 
z = 1) is perfectly consistent with the data. However, clusters 
make the confidence region substantially narrower (improve 
Wa constraints for a fixed Wq). A cosmological constant model 
(wo = -l,Wa = 0) is still consistent with the data. 

Finally, we note that in either case, the degeneracy between 
Wo and Wa is almost linear, Wa = A + Bwo- For such degen- 
eracies, constraints on constant w are equivalent to those for 
evolving w a t the pivot redshift, = (1 + Zp)^^ = 1 + 1/ 



& Jain|2()04| . From the slopes of degeneracies in Fig.[T2j we 
find Zp Ri 0.29 without clusters and Zp « 0.38 when cluster in- 
formation is included. Therefore, our combined constraints on 
constant w (§|8} can also be interpreted as those for evolving 
w at this pivot redshift. 

9.2. Wo in Non-Flat Universe 

The final case we consider is constant iv in a non-flat uni- 
verse. The cosmological grid in this case is (flj^, D.^, Wq, /i, ffs) 
with the requirement that the dark energy density is Ox = 
I-Om-Oj:. The hkelihood is marginalized over Om. h, as, and 
the constraints on fl^ and Wq are shown in Fig. 13 Including 
clusters does not noticeably improve the measurement of fit; 



-0.9 



f -1.0 



-1.1 



-1.2 



-1.3 




Fig. 13. — Equation of state constrains from WMAP, BAO, SN la, and clusters 
in the case of non-flat universe. We find wo = - 1.03 ±0.06 and flj. = -0.008 ± 
0.009 with all the data combined. 



by far the most significant contribution to the Clk constraint 
is from combination of WMAP and BAO data ( |Komatsu et al.| 
2008| l. However, clusters do substantially improve the equa- 
tion of state measurement: Wq = -1.03 ± 0.06 to be compared 
with ±0.085 without clusters. A flat ACDM model {Q.^ = 0, 
Wo = -1) is still consistent with the data within 68% CL. 



10. SUMMARY AND CONCLUSIONS 

We presented constraints on the cosmological parameters 
from a new measurement on the galaxy cluster mass function 
in the redshift range z = 0- 0.9. All major sources of informa- 
tion contained in the cluster mass function — its overall nor- 
malization and slope at z = 0, and evolution at high redshifts 
— are determined with our new data with a higher statistical 
accuracy and smaller systematic errors than before. This leads 
to much improved and more reliable constraints on the cos- 
mological parameters. 

From the normalization of the mass ftinction estimated at 
low redshifts, we derive the parameter degenerate with D.m- 
(T8(Om/0.25)° *^ = 0.813±0.013 (stat) ±0.024 (sys). The slope 
of the low-z mass function is a measure of Ciuh: = 
0.184 ± 0.037; combined with the HST prior on h, this is an 
independent measurement of CIm = 0.255 ± 0.043. The matter 
density can be independently measured with our cluster data 
using evolution of the temperature function, yielding consis- 
tent results, Om = 0.30 ± 0.05 in a flat ACDM model and 
0.34 ± 0.08 in a general cosmology. 

Evolution of the mass functions between z = and 0.5 (me- 
dian redshift for our high-z sample) constrains CIa = 0.83 ± 
0.15 in non-flat ACDM cosmology, or the dark energy equa- 
tion of state parameter, Wq = -1.14 ± 0.21, in a spatially flat 
Universe. Inclusion of the information provided by our clus- 
ter data also significantly improves the equation of state con- 
straints obtained from combination of multiple cosmological 
datasets. For example, by combining the 5-year WMAP, most 
recent supernovae measurements, and detection of baryonic 
acoustic oscillations in the SDSS with our cluster data, we ob- 
tain Wo = -0.991 ±0.045 (stat) ±0.040 (sys); both the statistical 
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and systematic errors in the combined constraint are a factor 
of 1 .5 - 2 smaller than those without clusters. Including cluster 
information also improves results for an evolving equation of 
state parameter and for constant w in a non-flat universe. A 
spatially flat ACDM model is within the 68% CL interval from 
the best fit in all cases that we tested. 

A good agreement between the geometric and growth of 
structure-based measurements of w in principle can be used to 
place Umits on modified gravity theories which att empt to ex- 
plain cosmic acceleration without dark energy (e.g., Wang et al. 
12007). When self- consistent models of non-linear collapse in 
such theories become available, it sould be straightforward to 
use our cluster data in such tests also. 

Comparison of the power spectrum normalization at z = 
obtained from clusters with the ampUtude of the CMB fluctu- 
ations is a sensitive measure of the mass of light neutrinos. We 
constrain 2 '"v < 0-33 eV at 95% CL, at the expense of slightly 
weakening the measurement of Wq obtained assuming that the 
neutrino masses are negligibly small. 

" http : //hea-ww . harvard . edu/400d/CCCP 



To facihtate the use of our cluster results in our cosmolog- 
ical studies, we provide at the project WWW site" machine 
readable tables of the likelihood function computed on several 
cosmological grids. 
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